Moving 10-years analysis of Phosphate at Northeast Atlantic Ocean for each season: - winter: January-March, - spring: April-June, - summer: July-September, - autumn: October-December Every year of the time dimension corresponds to the 10-year centred average of each season. Decades span : - from 1963-1972 until 2005-2014 (winter) - from 1963-1972 until 2005-2014 (spring) - from 1964-1973 until 2005-2014 (summer) - from 1964-1973 until 2005-2014 (autumn) Observational data span from 1962 to 2014. Depth range (IODE standard depths): -3000.0, -2500.0, -2000.0, -1750, -1500.0, -1400.0, -1300.0, -1200.0, -1100.0, -1000.0, -900.0, -800.0, -700.0, -600.0, -500.0, -400.0, -300.0, -250.0, -200.0, -150.0, -125.0, -100.0, -75.0, -50.0,-40.0, -30.0, -20.0, -10.0, -5.0, -0.0 Data Sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVA (Data-Interpolating Variational Analysis) tool. GEBCO 1min topography is used for the contouring preparation. Analyzed filed masked using relative error threshold 0.3 and 0.5 DIVA settings. Signal to noise ratio and correlation length were optimized and filtered vertically and a seasonally-averaged profile was used. Logarithmic transformation applied to the data prior to the analysis. Background field: the data mean value is subtracted from the data. Detrending of data: no, Advection constraint applied: no. Units: umol/l

The analysis was performed per season using DIVA software tool (Data-Interpolating Variational Analysis). The analyses products are stored as NetCDF CF files and made available as WMS layers for easy browsing and adding. Every step of the time dimension corresponds to a 6-year moving average from 1983 to 2016. The depth dimension spans from surface to 1000 m, with 21 vertical levels. The boundaries and overlapping zones between these regions were filtered to avoid any unrealistic spatial discontinuities. This combined water body dissolved oxygen concentration product is masked using the relative error threshold 0.5. Units: µmol/l Created by 'University of Liège, GeoHydrodynamics and Environment Research (ULiège-GHER)'. The data used as input for DIVA have been extracted from the EMODnet Chemistry Download Service: https://emodnet-chemistry.maris.nl/search Intermediate regional data products: Mediterranean Sea - DIVA 4D 6-year analysis of Water body phosphate 1968/2016 v2018, Arctic Ocean - DIVA 4D 6-year analysis of Water body phosphate 1980/2017 v2018, North Sea - DIVA 4D 6-year analysis of Water body phosphate 1980/2017 v2018, Black Sea - DIVA 4D 6-year analysis of Water body phosphate 1990/2016 v2018, North East Atlantic Ocean - DIVA 4D 6-year analysis of Water body phosphate 1960/2017 v2018, Baltic Sea - DIVA 4D 6-year analysis of Water body phosphate 1980/2016 v2018

Moving 10-years analysis of nitrate plus nitrite at Northeast Atlantic Ocean for each season: - winter: January-March, - spring: April-June, - summer: July-September, - autumn: October-December. Every year of the time dimension corresponds to the 10-year centred average of each season. Decades span : - from 1984-1993 until 2005-2014 (winter) - from 1979-1988 until 2005-2014 (spring) - from 1982-1991 until 2005-2014 (summer) - from 1972-1981 until 2005-2014 (autumn) Observational data span from 1962 to 2014. Depth range (IODE standard depths): -3000.0, -2500.0, -2000.0, -1750, -1500.0, -1400.0, -1300.0, -1200.0, -1100.0, -1000.0, -900.0, -800.0, -700.0, -600.0, -500.0, -400.0, -300.0, -250.0, -200.0, -150.0, -125.0, -100.0, -75.0, -50.0,-40.0, -30.0, -20.0, -10.0, -5.0, -0.0 Data Sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVA (Data-Interpolating Variational Analysis) tool. GEBCO 1min topography is used for the contouring preparation. Analyzed filed masked using relative error threshold 0.3 and 0.5 DIVA settings. Signal to noise ratio and correlation length were optimized and filtered vertically and a seasonally-averaged profile was used. Logarithmic transformation applied to the data prior to the analysis. Background field: the data mean value is subtracted from the data. Detrending of data: no, Advection consraint applied: no. Units: umol/l

Moving 6-year analysis and visualization of DIN in the North Sea. Four seasons (December-February, March-May, June-August, September-November). Data Sources: observational data from SeaDataNet/EMODnet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVAnd (Data-Interpolating Variational Analysis) tool, version 2.7.9. results were subjected to the minfield option in DIVAnd to avoid negative/underestimated values in the interpolated results; error threshold masks L1 (0.3) and L2 (0.5) are included as well as the unmasked field. The depth dimension allows visualizing the gridded field at various depths.

Moving 6-year analysis of Silicate at Atlantic Sea for each season. - winter: January-March, - spring: April-June, - summer: July-September, - autumn: October-December Every year of the time dimension corresponds to the 6-year centred average of each season. 6-year periods span - from 1972-1977 until 2015-2020 (winter), - from 1984-1989 until 2015-2020 (spring), - from 1972-1977 until 2015-2020 (summer), - from 1971-1976 until 2015-2020 (autumn). Observational data span from 1971 to 2020. Depth range (IODE standard depths): -2000.0, -1750, -1500.0, -1400.0, -1300.0, -1200.0, -1100.0, -1000.0, -900.0, -800.0, -700.0, -600.0, -500.0, -400.0, -300.0, -250.0, -200.0, -150.0, -125.0, -100.0, -75.0, -50.0,-40.0, -30.0, -20.0, -10.0, -5.0, -0.0 Data Sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVA (Data-Interpolating Variational Analysis) tool. GEBCO 1min topography is used for the contouring preparation. Analyzed filed masked using relative error threshold 0.3 and 0.5 DIVA settings. Correlation length was optimized and filtered vertically and a seasonally-averaged profile was used. Signal to noise ratio was fixed to 1. Logarithmic transformation applied to the data prior to the analysis. Background field: the data mean value is subtracted from the data. Detrending of data: no, Advection constraint applied: no. Units: umol/l

Moving 6-year analysis of Water body chlorophyll-a in the NorthEast Atlantic for each season: - winter: January-March, - spring: April-June, - summer: July-September, - autumn: October-December. Every year of the time dimension corresponds to the 6-year centred average of each season. 6-year periods span from 1971/1976 until 2016/2021. Observation data span from 1971 to 2021. High-frequency observation trails were filtered to a 3h temporal resolution. Depth levels (IODE standard depths): [0.0, 5.0, 10.0, 20.0, 30.0, 40.0, 50.0, 75.0, 100.0, 125.0, 150.0, 200.0, 250.0, 300.0]. Data sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Descrption of DIVAnd analysis: the computation was done with DIVAnd (Data-Interpolating Variational Analysis in n dimensions), version 2.7.4, using GEBCO 30 sec topography for the spatial connectivity of water masses. The horizontal resolution of the produced DIVAnd maps is 0.1 degrees. Horizontal correlation length varies from 200km in open sea regions to 50km at the coast. Vertical correlation length is defined as twice the vertical resolution. Signal-to-noise ratio was fixed to 1 for vertical profiles and 0.1 for time series to account for the redundancy in the time series observations. A logarithmic transformation (DIVAnd.Anam.loglin) was applied to the data prior to the analysis to avoid unrealistic negative values. Background field: a vertically-filtered profile of the seasonal data mean value (including all years) is substracted from the data. Detrending of data: no, advection constraint applied: no. Units: mg/m3.

Moving 10-years analysis of Silicate at Northeast Atlantic Ocean for each season. - winter: January-March, - spring: April-June, - summer: July-September, - autumn: October-December Every year of the time dimension corresponds to the 10-year centred average of each season. Decades span - from 1963-1972 until 2005-2014 (winter) - from 1979-1988 until 2005-2014 (spring) - from 1990-1999 until 2005-2014 (summer) - from 1964-1973 until 2005-2014 (autumn) Observational data span from 1962 to 2014. Depth range (IODE standard depths): -3000.0, -2500.0, -2000.0, -1750, -1500.0, -1400.0, -1300.0, -1200.0, -1100.0, -1000.0, -900.0, -800.0, -700.0, -600.0, -500.0, -400.0, -300.0, -250.0, -200.0, -150.0, -125.0, -100.0, -75.0, -50.0,-40.0, -30.0, -20.0, -10.0, -5.0, -0.0 Data Sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVA (Data-Interpolating Variational Analysis) tool. GEBCO 1min topography is used for the contouring preparation. Analyzed filed masked using relative error threshold 0.3 and 0.5 DIVA settings. Signal to noise ratio and correlation length were optimized and filtered vertically and a seasonally-averaged profile was used. Logarithmic transformation applied to the data prior to the analysis. Background field: the data mean value is subtracted from the data. Detrending of data: no, Advection consraint applied: no. Units: umol/l

Moving 6-year analysis and visualization of Water body chlorophyll-a in the North Sea. Four seasons (December-February, March-May, June-August, September-November). Data Sources: observational data from SeaDataNet/EMODnet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVAnd (Data-Interpolating Variational Analysis) tool, version 2.7.9. results were subjected to the minfield option in DIVAnd to avoid negative/underestimated values in the interpolated results; error threshold masks L1 (0.3) and L2 (0.5) are included as well as the unmasked field. The depth dimension allows visualizing the gridded field at various depths.

Moving 6-year analysis of Water body dissolved inorganic nitrogen in the NorthEast Atlantic for each season: - winter: January-March, - spring: April-June, - summer: July-September, - autumn: October-December. Every year of the time dimension corresponds to the 6-year centred average of each season. 6-year periods span from 1971/1976 until 2016/2021. Observation data span from 1971 to 2021. Depth levels (IODE standard depths): [0.0, 5.0, 10.0, 20.0, 30.0, 40.0, 50.0, 75.0, 100.0, 125.0, 150.0, 200.0, 250.0, 300.0, 400.0, 500.0, 600.0, 700.0, 800.0, 900.0, 1000.0, 1100.0, 1200.0, 1300.0, 1400.0, 1500.0, 1750.0, 2000.0]. Data sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Descrption of DIVAnd analysis: the computation was done with DIVAnd (Data-Interpolating Variational Analysis in n dimensions), version 2.7.4, using GEBCO 30 sec topography for the spatial connectivity of water masses. The horizontal resolution of the produced DIVAnd maps is 0.1 degrees. Horizontal correlation length varies from 400km in open sea regions to 50km at the coast. Vertical correlation length is defined as twice the vertical resolution. Signal-to-noise ratio was fixed to 1 for vertical profiles and 0.1 for time series to account for the redundancy in the time series observations. A logarithmic transformation (DIVAnd.Anam.loglin) was applied to the data prior to the analysis to avoid unrealistic negative values. Background field: a vertically-filtered profile of the seasonal data mean value (including all years) is substracted from the data. Detrending of data: no, advection constraint applied: no. Units: umol/l.

Combined product of Water body chlorophyll a based on DIVA 4D 10-year analysis on five regions : Northeast Atlantic Ocean, North Sea, Baltic Sea, Black Sea, Mediterranean Sea. The boundaries and overlapping zones between these five regions were filtered to avoid any unrealistic spatial discontinuities. This combined Water body chlorophyll a product is masked using the relative error threshold 0.5. Units: mg/m^3. Created by 'University of Liege, GeoHydrodynamics and Environment Research (ULg-GHER)'.